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Abstract 

We have studied the magnetic behavior of dextran-coated magnetite (Fe304) nanoparticles with 
median particle size (d) = 8 nm. Magnetization curves and in-field Mossbauer spectroscopy mea- 
surements showed that the magnetic moment Ms of the particles was much smaller than the bulk 
material. However, we found no evidence of magnetic irreversibility or non-saturating behavior 
at high fields, usually associated to spin canting. The values of magnetic anisotropy K e ff from 
different techniques indicate that surface or shape contributions are negligible. It is proposed that 
these particles have bulk-like ferrimagnetic structure with ordered A and B sublattices, but nearly 
compensated magnetic moments. The dependence of the blocking temperature with frequency and 
applied fields, Tb(H,uj), suggests that the observed non-monotonic behavior is governed by the 
strength of interparticle interactions. 
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* Present address: Instituto de Nanociencia de Aragon, Universidad de Zaragoza, Zaragoza, Spain 
^Electronic address: goya@unizar.es 



1 



I. INTRODUCTION 



Iron oxide nanoparticles have gained technological significance after they became known 
as media contrast agents in clinical protocols for Magnetic Resonance Imaging (MRI) and 
Magnetic Fluid Hyperthermia (MFH). For these applications, Magnetite (Fe3C>4 ) and 
maghemite (7-Fe20s) nanoparticles are the preferred phases because of their low toxicity, 
high saturation magnetic moment at room temperature (Ms ~ 75 and 82 emu/g, respec- 
tively) and the highest ordering temperatures among spinel ferrites.jfl] To improve stability 
and biocompatibility the particles are coated with a polysaccharide and suspended in water- 
based solvents. Despite the success attained in obtaining efficient dispersions of iron oxide 
particles for MRI protocols, the understanding of the effects of these coatings on the efficacy 
in MRI or MFH applications (from which a material with desired properties can be designed 
for specific purposes) is far from complete. 3, 0, U, 0] 

Regarding the fundamental mechanisms that govern the magnetic behavior of a magnetic 
colloid, the connection between single-particle properties and collective behavior of the en- 
semble of magnetic nanoparticles has many subtle facets. When entering the few-nanometer 
scale the magnetic behavior of a single particle gradually differs from the corresponding 
bulk materials, making increasingly inaccurate the description of the collective behavior in 
terms of 'bulk-like' entities. Dipolar magnetic interactions between ferromagnetic nanopar- 
ticles are known to affect the magnetic dynamics of the system through changes on the 
average anisotropy energy barriers, which in turn determines the relaxation times, magnetic 
hardness, and ordering temperature. Indeed, the effects of interparticle interactions are not 
restricted only to the anisotropy barriers but have also been proposed as stabilizers of the 
magnetic order at both particle core and surface p. In spite of sustained efforts to solve 
this problem alongthe last decade, the exact mechanisms linking these properties are still 
being discussed |7J,|8|. Theoretical approaches usually start from the original model proposed 
by Stoner and Wolfarth [9] for non-interacting, monodispersed single-domain particles, and 
add some specific perturbation to the collective or single-particle properties within different 
mathematical landscapes. In this way microscopic mechanisms like spin disorder, surface 



contributions and col 
original model. 



10 



ective behavior for strongly correlated particles have been added to the 
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121 ] On the other side, suitable experimental systems for testing those 



models are still a key problem to be solved, because of the difficulty of synthesizing samples 
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with different (controlled) particle size distribution and particle interactions. 

Magnetite is a mixed (Fe 2+ and Fe 3+ ) iron oxide with spinel structure, composed by 
a cubic close packed oxygen array, plus six (2 + 4) interpenetrating f.c.c. lattices of two 
nonequivalent cation sites with fourfold (8a) and sixfold (16d) oxygen coordination, la- 
beled as A and B sites, respectively. The cubic unit cell has lattice parameter a = 8.39 
A. [l(| Magnetic interactions within A and B sublattices are of ferromagnetic type, whereas 
the strong antiferromagnetic A-B coupling yields a high Curie temperature Tq = 850 K. 
Bulk Fe3C>4 is magnetically soft, i.a, the magnetization can be fully saturated in fields H < 
IkOe for any crystal direction. [12J At room temperature, the combination of first- and 
second-order magnetocrystalline anisotropy constants K\ = — 13 x 10 4 erg/cm 3 and K 2 = 
—3 x 10 4 erg/ cm 3 makes the easy axis of magnetization to be along the crystallographic [111] 
direction. [18( At low temperatures magnetite undergoes a first-order transition, namely the 
Verwey transition, located at Ty ~ 115 — 124 K depending on sample conditions. [19] The 
most conspicuous evidence of this transition is the sharp increase of the resistivity by two 
orders of magnitude, 19., 0] but also other properties such as magnetization, thermopower 
and reflectivity show abrupt changes. Regarding the structural transition reported at Ty, 
there is some consensus about the ortorhombic symmetry of the low-T phase, but no con- 



clusive structural model has yet emerged. [2l|, |22j Despite the different crystal structures 
proposed for the low temperature phase, magnetic measurements have shown that below Ty 
the magnetocrystalline anisotropy energy can be described as having orthorhombic symme- 
try with the easy c axis tilted ~ 0.2° from the <1 0> cubic axis, and the new a (hard) and 
b axes lying on the former cubic (001) plane, The resulting anisotropy is then uniaxial 
with constants K a « 2 x 10 6 erg/cm 3 and 21 x 10 4 erg/cm 3 , considerably larger than 



the K\ constant of the high-T phase. 
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In this work we performed a detailed magnetic and structural characterization of a col- 
loidal suspension of dextran-coated magnetite nanoparticles, with the purpose of systemat- 
ically characterize both the single particle properties and the resulting collective behavior. 
The present paper is organized as follows: in Sec. |Hj the experimental procedures are 
described. In Sec. IIHI the results are presented and discussed in four parts: Sec. IIII Al 
deals with the analysis of the structural data. In Sec. IIII Bl the results from magnetization 
measurements are shown. Sec. IIII CI describes the Mossbauer data obtained with different 
applied fields and temperatures. Finally, the discussion and conclusions drawn from the 
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present work are given in Sec. IIV1 



II. EXPERIMENTAL PROCEDURE 

The samples studied in this work consisted of Dextran-coated magnetite nanoparticles 
dispersed in a water-based ferrofluid (ENDOREM™, Guerbet) used in clinical protocols as 
a contrast agent for magnetic resonance imaging of the liver and the spleen. j^J The disper- 
sion of nanoparticles was studied as supplied, i.e. suspended in an aqueous liquid carrier, 
conditioned in sealed cylindrical sample holders (5 mm diameter x 6 mm high). A second 
fraction of the sample was lyophilized (@ T = 85 K and P = 10~ 4 Torr) to increase the 
particle concentration keeping the particle size distribution unaltered. These as supplied and 
lyophilized samples will be labeled hereafter as END1 and ENDS, respectively. Characteri- 
zation through X-ray diffraction (XRD) was performed in a Philips PW 1820 diffractometer 
using Cu— K a radiation (A = 1.5418 A) and Ni filter. Transmission Electron Microscopy 
(TEM) images were performed using a 200 kV Model CM200 Philips electron microscope, 
conditioning the samples by dropping an alcohol-powder suspension on a carbon-coated 
nickel grid. The elemental analysis of Carbon, Hydrogen, Nitrogen and Sulphur contents 
was performed with a Perkin-Elmer 2400 CHNS microanalizer. The Mossbauer spectroscopy 
(MS) measurements were performed between 4.2 and 296 K in a liquid He flow cryostat, 
with a conventional constant-acceleration spectrometer in transmission geometry using a 
57 Co/Rh source with c.a. 50 mCi activity. For in-field Mossbauer measurements, the pow- 
der samples were prepared between acrylic discs, and mounted in the bore of a 140-kOe 
superconducting magnet, in a vertical source-sample-detector setup such that the direction 
of gamma ray propagation was parallel to the magnetic-field axis. For this setup, a sine- 
shaped velocity waveform was used to minimize mechanical noise. The spectra were fitted 
to Lorentzian line shapes using a non-linear least-squares program, calibrating the velocity 
scale with a foil of a-Fe at 296 K. When necessary, a distribution of hyperfine magnetic 
fields, isomer shift and quadrupole splitting have been used to fit the spectra. Magnetiza- 
tion and ac magnetic susceptibility measurements were performed in a commercial SQUID 
magnetometer both in zero-field-cooling (ZFC) and field-cooling (FC) modes, between 1.8 
K < T < 250 K and under applied fields up to H = 70 kOe. The frequency dependence of 
both in-phase x'(T") and out-of-phase x"(T") components of the ac magnetic susceptibility 
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were measured by using excitation fields of 1 — 4 Oe and driving frequencies 0.01 Hz < / < 
1500 Hz. 



III. EXPERIMENTAL RESULTS 



A. Structural Analysis 



The XRD profile of the lyophilized powder (fig. 1) was composed of broad lines that could 
be indexed with a cubic spinel crystal structure. The presence of minor amounts of unknown 
organic ingredients (used for coating and stabilizing the ferrofluid) were also detected. A 
rough estimation of the average grain size (d) RX = 7 ±1 nm was obtained by applying the 
Scherrer formula to the most intense XRD line of the spinel structure, without including 
effects from crystal stress. This is in agreement with previously reported low-temperature 



Mossbauer data [221 I2S( that showed the existence of Fe3C>4 as the only phase, without 
evidences of the related 7-Fe2C>3 (maghemite) or a-Fe2C>3 (hematite) phases. 

Fig. 1 also shows the particle size histogram obtained by measuring the diameter of about 
200 particles in Transmission Electron Microscopy (TEM) images. The resulting histogram 
could be fitted with a log-normal distribution, with median diameter (d) TEM = 8.3 ± 0.8 
nm and distribution width = 0.6. It can be noticed from comparison between the fitted 
curve and the experimental histogram that the counted number of particles with d < 3 nm is 
considerably larger than the expected from a pure log-normal size distribution with the fitted 
parameters. The magnified TEM image of specific areas showed particles with spherical-like 
shapes, although for few particles with the largest diameters we observed faceted particles, as 
selected in the inset of Fig. 1. The latter shape is probably related to the high crystallinity 
of the particles reflecting the cubic crystal habit of the Fe 3 4 phase. 

B. Magnetization Data 

Some of the basic magnetic features of the present particles established previously [3, 13| 
are summarized here for clarity. j^J The particles were superparamagnetic (SPM) at room 
temperature, and magnetization curves M(T) taken in FC and ZFC modes showed blocking 



temperatures T B of ca. 50if (for H = 100 Oe). [2jj Additionally, ac susceptibility data 
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x{T, f) confirmed |28( the thermally activated (Arrhenius) nature of the blocking process, 
and indicated that the effect of interparticle interactions was to increase the energy barriers. 

In order to measure the evolution of the saturation magnetization of the particles with 
temperature, we applied a field of H = 70 kOe at T = 300 K (i.e., above the melting 
temperature of the liquid carrier) to the as provided sample, cooled down the system and 
measured the Mpc{T) curve for increasing temperatures keeping the same field. It can be 
seen from fig. 2 that these fully aligned particles show a nearly linear dependence of the 
Mpc{T) data, varying within ^13 — 7 emu/g in the 1.8 < T < 300 K temperature range. 
All magnetization values were converted into emu/g Fes04 units by subtracting the mass 
of organic components as determined from elemental CHNS analysis, calibrated within 0.1 
% precision (Table HJ. These results were in agreement with the expected nanoparticle con- 
centration from the nominal composition of the ferrofluid (0.07 mol Fe3C>4 /l) and yielded 
saturation values Ms = 10-12 emu/g, much smaller than the expected 85-95 emu/g for 
bulk magnetite. Additionally, the M(T) curves did not show any indication of the jump in 
the magnetization expected for the Verwey transition. |19J] To explore possible aging (oxida- 
tion) effects we compared our initial measurements with other runs performed after several 
months, taken in three different magnetometers, finding that the magnetic moment was re- 
produced within few percent of precision. Additional support for the observed reduction of 
the magnetic moment of the particles came from the M(H) curves both below and above the 
blocked temperature Tb ~ 40 K (see below). As will be discussed in Sec. IIII C| Mossbauer 
data showed no indication of other iron-oxides of smaller magnetic moment, as could be ex- 
pected in case of particle oxidation (hematite Fe20s) or hydration (ferrihydrite). Therefore 
we attempted to verify these rather low values by independent measurements, namely, the 
magnetic behavior in the SPM state. 

The magnetization Mspm{H,T) of an ensemble of non-interacting, single-domain parti- 
cles in the SPM state, each with magnetic moment /i, can be represented by a Langevin 
function L(x) of argument x = fiH/ksT. Therefore all data M(H,T) should scale into a 
single curve when plotted against H/T, provided that Tb < T C Tg. Based on the log- 
normal distribution of particle sizes obtained from TEM d ata, we modify the model above 



by using a distribution-weighted sum of Langevin functions 



30] 
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where f(fi) is the log- normal distribution function 
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of magnetic moments \x. In this equation a M is the distribution width and (//) is the median 
of the distribution related to the mean magnetic moment fi m by \i m = (u) exp(cr^/2). j^l 

The best fit for sample END1 (for which dipolar interactions are expected to be less 
important) yielded a median magnetic moment (//) = 290 fiB and a distribution width = 
1.68 for the highest measured temperature (see figure 3). The fittings of M(H) curves were 
performed for T = 100, 150, 200 and 250 K. We have observed that the resulting median value 
(/j) increased c.a 10% from 250 K to 100 K, which is likely to be related to the temperature 
dependence of the magnetic moment M(T) discussed above (fig. 2). The distribution width 
(Tp decreased ~ 6% within the same temperature range. To relate this moment distribution to 
a particle size distribution we have assumed spherical geometry for the particles, as observed 
from TEM data (i.e., (//) = Mg| (d) 3 ) and the experimental saturation magnetization value 
Ms = 7.2 emu/g from the M(H) data at the same temperature. We have obtained (d) MAG = 
7.7 ± 0.5 nm with <Jd = 0.56. The median particle size {d) MAG obtained with the measured 
Ms value is in close agreement with the (d) TEM = 8.3 ± 0.8 nm from TEM micrographs, 
and the corresponding distribution widths are identical within experimental accuracy. The 
outsized reduction in the magnetic moment of the particles observed in M(T) curves (fig. 2) 
in the blocked state is therefore supported, independently, by the measured M(H, T) curves 
in the SPM state. The origin of this reduction will be discussed in Sec. IIVI We note here 
that the latter procedure for obtaining the particle magnetic moment fx does not depend on 
the absolute magnetization values but on the fiH product for a given temperature. 

We have further characterized the transition from blocked to SPM state by measuring the 
hysteresis curves M(H) at different temperatures between 1.8 and 100 K for samples END1 
and ENDS. From the obtained M(H, T) curves (fig. 4) it is observed that the coercivity fields 
He of both samples remain essentially zero down to T 17 — 18 K, increasing appreciably 
below this temperature. The temperature dependence of H C (T), as extracted from M(H, T) 
curves, is the expected for single-domain particles, i.e., 
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He = H C o 



1 - 



2hk B T\ 3 

KeffV 



(3) 



where Hco = 2aK e ff /Ms, and a = 0.48 is a phenomenological constant used for a randomly 
oriented distribution of easy axis. [^] 

Therefore, eqE] can be used for determining K e ff from the slope of the Hc(T) curves, 
through the product E a = K e ffV, using the extrapolated value H C q. Expressing E a in terms 
of an anisotropy field E a = H^Ms, the effective anisotropy constant K e ff for a cubic system 
with easy axis along the [111] direction is related to K\ and K 2 through K e ff = + ^K 2 , 
[31I ] and thus K e ff = 18.7 x 10 4 erg/cm 3 for bulk Fe3C>4 . The values of K e ff were extracted 
from the best fit of eq. El for T < 15 K, and the calculation was made using both the mean 
value from TEM distribution (d) = 8.3 nm and the experimental Ms = 12 emu/g values 
from M(H) data at low temperatures. The resulting K e ff values shown in Table HP are very 
close to the expected anisotropy of bulk magnetite, suggesting that no major contributions 



so in agreement with our 
28j |. By extrapolating the 



other than magnetocrystalline are important. These values are a 
previous results obtained from low temperature Mossbauer data 
linear fits to Hc{T) = 0, we found values of Tb = 17 A K and 18.5 K for samples END1 and 
ENDS, respectively. We mention that these values obtained from the extrapolation of Hc(T) 
are smaller than the blocking temperatures obtained from ZFC curves, which we associate to 
the different field regimes of each measurement (see discussion in Sec. IIV|) . The extrapolated 
T — > K values were Hqo = 418 Oe and 428 Oe for samples END1 and ENDS, respectively. 
It is worth noting from fig. 4 that the saturation of the magnetization (M > 0.95 Ms) is 
attained for H w 10 kOe, the value of Ms being calculated from the H~ l — > extrapolation 
in the M vs. H~ l curves. These saturating fields are low when compared with results 
from other nanoparticles of similar avera ge s izes, where non-saturating behavior extended 



for applied fields larger than 70-90 kOe.[32, |33j, |34j For sample END1 the extrapolated 



Ms = 13.5 emu/g value from the M(H) curves at T = 1.8 K (fig. 4) is slightly larger than 
the Ms = 12.9 emu/g value from FC data (fig. 2). The origin of this difference is probably 
related to differences in particle alignment when the field is applied in the liquid state (FC 
process) and below the blocking temperature (M(H) curves). 
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C. In-Field Mossbauer Data 



The reversal of the magnetic moment of single-domain particles is governed by the rate 
of flipping over an energy barrier that is proportional to the volume of the particle. This 
thermally activated process is described by a relaxation time r that follows an Arrhenius 
law 



r=r ° exp GSf)' (4) 

where E a is the energy barrier for the magnetic reversion of the system, fcg is the Boltzmann 
constant, and r is the pre- exponential factor (of the order of 10"' - 10"" sJ.QQ At room 
temperature (T ^> Tg) the relationship r << tl (where tl = 10~ 8 — 10~ 9 s is the Larmor 
precession time of the nuclear magnetic momentsj^]) is verified, so the fast relaxation 
makes the magnetic interactions to be averaged to zero and thus yields a paramagnetic-like 
doublet. For temperatures below Tg (i.e., r >> tl), the Mossbauer spectrum is magnetically 
split. The hyperfine parameters at low temperature reported in ref. were in agreement 



38, 



391 and showed no sizeable reduction of 



with other reports on nanostructured magnetite 

' a n i 

hyperfine fields, in contrast to other nanoparticle systems |4d.l4lLl42l l43| where the reduction 
effect was assigned to surface disorder effects. From eq. 0] it follows that the blocking 
temperature Tg depends on the experimental window time tm of a given measurement. It 
is well known that the ratio between the blocking temperatures from different techniques 
like Mossbauer and magnetization measurements (having typically tm = 10 -8 and 10 2 s, 

^ rpmOSS 

respectively) inferred from eq. El can be as large as JL as ~ xo. The blocking temperature 

1 B 

for Mossbauer experiments are usually defined as the temperature at which each of the 
blocked and SPM subspectra have 50% of the total spectral area. We estimated Tg oss from 
a series of Mossbauer spectra taken for 4.2 K < T < 294 K (not shown) as being located 
somewhere between 110 and 120 K, which is about 2 — 3 times the value obtained from 
magnetization data. 

In order to obtain local information of the rotational process of the magnetic moments, 
we performed in-field Mossbauer measurements up to B app = 120 kG. In this kind of exper- 
iments the effective field B e ff actually observed is the vectorial sum of B app and the internal 
hyperfine fields Bh yp , B e ff = Bh yp + B app thus providing useful information about relative 
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orientation between the local magnetic field at the probe nucleus and the external field (i.e., 
between B app and the magnetic sublattices). Figure 5 shows the resulting Mossbauer spectra 
of sample ENDS at 4.2 K for B app = 0, 10, 30, 50, 80, and 120 kG, applied along the 7-ray di- 
rection. It can be seen from the figure that for increasing applied fields two distinct spectral 
components are apparently resolved. These two components originate from the ferrimag- 
netic alignment of Fe atoms at A and B sites in the spinel structure, and the assignment of 
each subspectra to A and B sites (see Table IIII|) is based on the zero-field spectrum. 

It can be observed in figure 5 that the intensities of the second and fifth lines change 
as the applied field increases, which is a direct indication of the progressive reversal of the 
spins in the direction of the field. The spectral intensities of each line of a magnetically split 
spectrum are usually expressed as 3 : x : 1 : 1 : x : 3, where x is the intensity of lines 2 or 
5. For the present setup, the values of x span from x = 2 for a randomly oriented system 
to x = for full alignment in the direction of the external field. More generally, the ratio of 
spectral areas (Ai) of lines 1 and 2 (or 5 and 6), can be related to the angle (3 between the 
magnetic quantization axis (assuming that the magnetic interaction is dominant) and the 
direction of the 7-ray beam (identical to the applied field B app in our experimental setup) 
through the expression 



sin 2^ = _^_ ; (5) 
H 4 + 3A' v ; 

where A = A%jA\ (or A^/Aq). The resulting values of A and j3 obtained from the fits at 
different fields are displayed in Table UTTl For B app = both subespectra could be fitted with 
the A = 0.67 value expected for a random powder. On the other side, we observed different 
behavior for A and B sublattices under increasing applied fields. Indeed, Table 11111 shows 
that the magnetic moments at A sites are already aligned for B app > 10 kG, whereas the 
corresponding B-site moments show a gradual rotation and remain slightly misaligned up to 
the largest field of B app = 120 kG. It is worthwhile to mention that the above observation 
relies on the capability of separately fitting (i.e., experimentally resolve) lines 2 and 5 of each 
subspectrum. These lines usually display strong overlap in ferrimagnetic materials, yielding 
unreliable A values from the fitting procedure. For the present data, we were able to resolve 
clearly lines 2 and 5 from each sublattice at B app = 120 kG, using the obtained values of 
intensities as the starting free parameters for fitting at lower B app values. 
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It is known that ferrimagnetic nanoparticles display non-collinear (canted) spin structure, 
although the location o 



being discussed. 



44 



45 



these disordered spins (i.e., at the surface, core, or both) is still 
4|| The actual spin structure is known to be more complex than 
the expected from the original Yafet-Kittel (YK) model of competing interactions, 47 1 
and neutron diffraction and Mossbauer studies have provided evidence that spin canting 



can be either restricted to a single magnetic sublattice or extended to both cation sites. 
3, 0,0, Eii- The present in-field Mossbauer data clearly show that magnetic moments 
at A-sites (and also a fraction at B-sites) are magnetically 'soft', in agreement with the 
saturating behavior observed in M(H) (see fig. 4) at low temperatures. Additionally, the 
resolution of the two magnetic sublattices due to the application of a large external field 
makes possible to identify a minority of magnetically hard spins as located at B-sites in the 
spinel structure. Despite their very high local anisotropy, however, this small fraction does 
not seem to contribute appreciably to the global magnetic anisotropy, as observed from the 
magnetization data presented before. 



IV. DISCUSSION 



The reduction of magnetic moments in small particles (1 — 10 nm) and thin films 
has been well documented along the last years 36:, 52] in both metallic and insulating 
materials 00000. A very-well accepted explanation for the reduction of Ms 
in ultrafine maghemite particles was given by Berkowitz et al. 53], who proposed a shell-core 
structure with a magnetically dead surface layer originated from demagnetization of the 
(paramagnetic) surface spins. If the above mechanism is assumed for the present particles, 
a simple calculation using the bulk magnetization value (M s ulk = 90 emu/g) would imply 
a core diameter (d) core = 2.5 nm and shell thickness t = 2.9 nm. Such a thick paramag- 
netic surface layer, although not strictly unphysical, appears very unlikely since Mossbauer 
spectra ruled out secondary Fe-containing phases, and CHNS compositional data showed 
that the organic material is much less than the necessary for such an organic coating. More 
recently, spin canting has been proposed ^,54, 0] as the mechanism for the Ms reduc- 
tion in spinel nanoparticles, due to competing interactions between sublattices that yield 



magnetic disorder at the particle surface. 



.56] 



Experimental evidence for this model 



found from Mossbauer spectroscopy, and neutron diffraction measurements. 
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tremely low values of Ms have been measured in ferrimagnetic spinel particles with d < 15 
nm, attaining values as small as ~ 0.06M| u ' fc , where M s ulk are the corresponding values of 



the bulk material. 



56 



58, 



59|. From the theoretical side, Monte Carlo simulations using 
IsoL Ig1 | have shown that the reduction of Mg 



different models and approximations |o^. I55L 
is size-dependent, and originates in a highly disordered spin surface with large magnetic 
anisotropy. The large magnetic anisotropy of the surface could be, in turn, responsible for 
the non-saturating M(H) curves even at large (H > 100 kOe ) applie d fields, in agreement 
with the observed behavior of many spinel nanoparticles. 43, 58, 62, 3| 

For present Fe3C>4 particles, however, the reduced Ms observed is not accompanied by 
the linear increase in M(H) at high fields characteristic of canted systems. On the contrary, 
both the M(H) curves and in-field Mossbauer spectra at low temperature showed that 
the magnetic saturation is attained at moderate (H ^10 kOe) fields. We are therefore 
induced to hypothesize that, instead a thick spin-disordered surface layer, the origin for 
the observed reduction in Ms could be related to internal compensation of the magnetic 
sublattices A and B in the spinel structure. Further support to this idea comes from the 
evolution of the relative intensities of the A and B subespectra for increasing fields: it can 
be seen from Table IIHI that the intensity of the subespectrum assigned to B site increases 
at expenses of the A subespectrum. The mechanisms that originate this peculiar situation 
are not obvious, and further experiments should be done before an explanation can be given. 



Turning now to the measured magnetic anisotropy of the present nanoparticles, we want 
to discuss the values extracted from different techniques involving different experimental 
conditions. It can be seen from Table |H] that the resulting values of E a obtained from 
magnetization and Mossbauer data are essentially coincident with the magnetocrystalline 
anisotropy of bulk Fe 3 4 . Indeed, this is a remarkable result since shape anisotropy is usually 
a major contribution to the average total anisotropy j^4|. It is generally accepted that, 
when the size of a magnetic particle decreases, surface effects become increasingly important. 
However, symmetry arguments suggest tha t the contribution from the surface should average 



to zero for a perfectly spherical particle. |40J] Moreover, numerical calculations have shown 
that for an unstrained simple cubic lattice the bulk anisotropy vanishes in first-order, whereas 
the contribution from the Neel surface anisotropy is second-order, having the same cubic 
symmetry of the bulk. Therefore for highly symmetric particle shapes such as cubic or 
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spherical, the contribution from both core and surface anisotropies should be comparable to 
the bulk magnetocrystalline values. These symmetry arguments are clearly applicable to our 
samples since, as mentioned above, the low magnetic anisotropy values obtained for these 
particles indicate magnetocrystalline anisotropy as the major contribution, setting an upper 
limit (no more than ~ 2 %) for possible deviations from sphericity. [64] Similar arguments 
can be applied to the case of cubic-shaped particles like the faceted ones observed in our 
TEM images, since demagnetizing factors for a rectangular prism are within 15 % of those 
of an ellipsoid having the same a=b, c, and r=c/a parameters, as calculated by Aharoni-Jf^l 
Therefore, the above results support a picture of particles with spherical or edge-rounded 
cubic shapes, having cubic crystal structure with bulk-like anisotropy. 

A third value for E a jk^ was previously j3| estimated from low temperature Mossbauer 



data, using the model of collective magnetic excitations (CME). |67( In this model, small- 
angle fluctuations of the magnetic moments around the easy magnetization axis, with char- 
acteristic times much shorter than the Mossbauer sensing time 7 7,, y ield a decrease of the 
observed magnetic hyperfme field B e ff(T) of Mossbauer spectra. |3j| Expressing the mag- 
netic energy E = E(a x , a y , a z ) of a single particle as a function of the direction cosines 
(a x , a y ,a z ) of the saturation magnetization Ms relative to the crystallographic axes, the 
thermal dependence of -B e // derived for an arbitrary function of the energy can be expressed 
asH 

BeffiT) = Bo[l~^fy (6) 
where Bq is the hyperfme field at T = 0, and the effective anisotropy constant k is given by 



1 1 

k~ ~ 2 



c^V 1 fd 2 E 



()( I ~ / V dOty 



(7) 



In a previous work on Fe3C>4 particles 



381 ] with median sizes (d) ranging from 150 to 5 nm, 
it was reported that for (d) < 50 nm no evidence of the Verwey transition is observed, sug- 
gesting that those particles preserve the cubic phase at low temperatures (T < Ty). There- 
fore, for the present nanometric Fes04 particles the values of magnetocrystalline anisotropy 
constants K\ and K 2 of the high-temperature cubic phase with easy magnetization axis [111] 
should be used at any temperature. Applying these considerations and eq. [7| an effective 
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anisotropy constant k = —^{K\ + ^j 2 -) is obtained. Using the (d) TEM value from the log- 
normal fit to the TEM distribution, the best fit of B e ff(T) gave an effective anisotropy k = 
24.3 x 10 4 erg/cm 3 , in good agreement with the 18.7 x 10 4 erg/cm 3 expected using the bulk 
magnetocrystalline values K\ and K 2 . 

Inasmuch as the collective magnetic excitation model refers to the thermal excitations 
of the lowest energy levels, for T << Tg, the E a parameter obtained from eq El using 
the CME model at zero field, gives information about single-particle anisotropy properties 
since interparticle interactions will have negligible effects on the detailed shape of the 
E a (9, H) surface in the vicinity and at the points of minima. On the other side, the values 
extracted from the frequency- dependent maxima in a.c. susceptibility data are related to 
the overcome of the magnetic anisotropy when thermal energy equals the average height 
of the energy barriers. The resulting E a values extracted from the unblocking process can 
therefore be influenced by other contributions (e.g. interparticle interactions) that affect 
the values of the E a (8, H) surface near the maxima. 62] 



In order to experimentally check the influence of small applied fields to the energy barri- 
ers, we measured the field- dependence of T B through ac susceptibility curves x{f-> T, H) at 
different applied dc fields. From the x{fi T, H) data taken with H = 10 and 150 Oe, we 
extracted the Ts{f,H) dependences shown in fig. 6. It can be observed that the frequency 
dependence T^ln/) is linear for both applied fields, demonstrating that the thermal acti- 
vation mechanism of the blocking transition is preserved. The best fits of T^" 1 (ln/) data 
were obtained with E a /ks = 436(9) and 296(15) K for H = 10 and 150 Oe, respectively, 
indicating that the effect of the applied field is to decrease the effective anisotropy barrier 
E a (H). These results give support to the idea that the application of small fields can perturb 
the relative values and position of the maxima of the anisotropy energy function. 

To further verify the above statements on a quantitative basis, we calculated the energy 
curves for a single particle in applied fields H = 0, 10 and 150 Oe by the following relation 



E(9, H) = K eff V + fx HM cos 6, (8) 

where 9 is the angle between the easy magnetization axis and the magnetization vector M 
of the particle, and H is the external field, applied along the easy axis of magnetization. In 
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the presence of a small applied field H parallel to the easy direction of the particle (assumed 
here to be uniaxial), the energies at the minima for 9 = 0° and 9 = 180° will shift, and the 
transition probabilities between the two minima will result in different relaxation times r + 
and r_ (i.e., for forward and backwards relaxation processes, respectively). In the present 
case of modest applied fields, the observed relaxation times can be expressed as 3^1 



111 

- = - + -, 9 

r r_ r + 

with 

r± = T exp(E(±H)/k B T). (10) 

The energy curve for H = 150 Oe was computed using a gradient method to calculate the 
local energy minima for each spin configuration along the magnetic energy surface E(9, H). 
Using eqs. IH1 - ITU1 and considering the experimental values M s = 12 emu/g, (d) TEM = 
8.3 nm and E a = K e ffV = 7.26 x 10 -14 erg the resulting height for the energy barrier 
AE/ks = Eg = go/kB — Eg = o/kB was 500 K at T = 10 K, in excellent agreement with the 
experimental results in Table ITT1 A systematic series of computation for increasing applied 
fields H from to 300 Oe showed that the difference AEq between the energy maximum and 
the minimum at 9 = goes deeper whereas the corresponding AE n (minimum at 9 = n) gets 
shallower. Therefore, the corresponding relaxation times r + and r_ show opposite behavior 
for increasing fields H provided that fioHM « K e ffV. We also computed the ratio t+/t_ 
for our experimental conditions. The results obtained for T = 10 K and the experimental 
applied fields H — 10 and 150 Oe (see fig. 6) yielded r + /r_ = 2 x lO^ 1 and 7 x 10~ 5 , 
respectively. These results give quantitative support to the assumption that the energy 
barriers (i.e., the whole relaxation process) can be changed in the observed amounts by the 
small fields used in our measurements. It should be pointed out that an exact calculation of 
the relaxation time would require, within Brown's model, to work with a probability function 
of angles and time, which obeys a Fokker-Planck equation. For symmetries different than 
uniaxial, the corresponding eigenvalues of the Fokker-Planck equation have been obtained 3(J 
by numerical methods only for asymptotical situations (i.e., ^7 — > 00 or — > 0) not 
applicable to our present case. 

It was already mentioned that the values of T# obtained from the extrapolation of 
Hq{T) — > from M(H) curves at low temperature differ from the ones obtained from 
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ZFC measurements. The difference between these two results is probably related to the 
different field regimes in both techniques: while the ZFC maxima are obtained at low fields 
(typically 100 Oe), from a previously demagnetized state, the Hq(Tb) is obtained after mag- 
netic saturation of the samples (we used H max = 30 kOe). Therefore, the moderate fields He 
applied when magnetization reversal takes place are likely to assist the unblocking process 
by decreasing the energy barriers and the inferred Tb will be lower than the correspond- 
ing values extracted from the ZFC curves measured at lower fields. The latter values, in 
turn, are mainly determined by the dynamics of the unblocking process that depends on 
the local minima at the energy surface originated from interparticle interactions. The above 
interpretation is in agreement with our estimations of the energy barriers E a /kB from a.c. 
susceptibility measured at low field, 2^1 (reproduced in Table ITT1 for clarity) which are larger 
than those extracted from the Hc(T) fit. 

During the analysis of the present data, the activation energies values extracted from 
x(f, T, H) data with H = and 150 Oe were compared to the previously reported value 
of E a = 694 K for the same END1 sample, obtained with H = 30 Oe (see Table ITTjl . 
Contrary to the expected, the value obtained at an intermediate dc field is larger than 
the two obtained at lower and higher fields. Having verified that the temperatures of the 
maxima in our x{f\ T) curves were reproducible within AT < 3 K on different experimental 
runs, we attempted to investigate this seemingly non-monotonic behavior of the anisotropy 
energy barriers by measuring the dependence Tb(H) from ZFC curves at different fields. 
The results are shown in figure 7. It can be seen that for sample END1 there is a continuous 
increasing of Tb(H) for decreasing fields, but a clear bending of Tb towards smaller values 
is observed for H < 20 Oe. The reproducibility of these maxima in M(T, H) curves was 
checked by performing a second set of ZFC curves, after allowing the sample to demagnetize 
at room temperature for several weeks. The results of the second run, plotted in fig. 7 (black 
dots) demonstrate that the main features of the Tb{H) curve obtained on the second run 
remain unaltered. This 'reentrant' behavior has been originally found in a rather different 
context, i.e., in CuMn and GdAl spin glasses, and interpreted as due to the existence of a 
finite transition temperature at zero field. [69] Theoretical evidence of reentrant behavior of 
the Gabay-Tolouse line has been also found using a classical m— vector model with uniaxial 

n 

anisotropy. [70] 

It seems more plausible, however, to relate this non-monotonic field dependence of Tb{H) 
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at low fields with changes in the anisotropy energy surface E(8,H). It can be seen from 
fig. 7 that the Tb{H) curves tend to merge each other in the large (> 500 Oe) region, in 
agreement with the erasing of local minima by the external field, and both curves steadily 
diverge for small applied fields. However, for the concentrated sample a monotonic behav- 
ior persists down to the smallest applied fields. This behavior has been already found in 
diluted magnetite-based ferrofluids with different concentrations, jzi], but no explanation 

'or the non- monotonic Tb{H) dependence was offered. In a detailed work, Sappey et al. 

~i\ proposed that the applied field can produce a small broadening of the energy barrier 
distribution, which suffices to give the observed field dependence of Tb- For the present 
data, different behavior is observed in two samples with identical particle size distributions 
but different average particle distances. A rough estimation using the compositional data 
from Table |U and (d) TEM = 8.3 nm, yielded values of the average particle-particle distance 
(e) = 63 nm and (e) = 2 nm for END1 and ENDS, respectively. Therefore in sample END1 
we have (e) ~ 8 (d) and dipolar interactions are likely to be much weaker than for the (e) ~ 
^j- situation of concentrated ENDS sample. The above arguments support the idea that the 
observed differences in Tb{H) curves are related to the influence of interparticle interactions 
on the energy barrier. The absence of nonmonotonic behavior in the concentrated sample 
indicates that the local dipolar field must be included when considering the influence of the 

n 

applied field on the barrier distribution, as proposed by Sappey et al. 72]. Whatever the 
actual mechanism of the observed behavior of Tb{H), the values of E a in Table HT1 obtained 
at H = 10, 30 and 150 Oe can be understood as having the same origin. 

In conclusion, our present results on nanometer-sized magnetite nanoparticles showed 
some intriguing facets. Single-particle anisotropy energy both in blocked and SPM regimes 
and different experimental time windows, indicates a reduced magnetization with minimum 
spin disorder. The magnetic response of the system can be described as 'bulk-like' 
(regarding magnetic softness and anisotropy) and, together with the reduced M s and 
the magnetic saturation in moderate fields experimentally observed, can be explained if 
a change in each sublattice magnetization is assumed. This in turn might be originated 
from departures from the bulk populations of the Fe 2+ and Fe 3+ ions at A and B sites 
yielding a more complex magnetic configuration than the two-sublattice structure. At low 
fields, particle interactions become increasingly relevant, and may erase the non-monotonic 
dependence of the blocking process observed in diluted samples. As the energy barrier 
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height and distribution usually determine the dynamics of magnetic nanoparticle systems 
at low fields, its complex dependence with external parameters has to be solved before 
accurate parameters can be extracted from magnetic measurements. 
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TABLE I: Sample composition obtained from CHNS elemental analysis of Carbon, Hydrogen, 
Nitrogen and Sulphur contents. The [NP] column is the calculated particle concentration (in 
particle/cm 3 ). 



Sample H (%) C (%) N (%) S (%) [NP] 

END1 60.3(5) 26.4(5) 3.3(5) 1.1(5) 6.5 x 10 

ENDS 2.0(5) 24.2(5) 0.8(5) 0.0(5) 1.6 x 10 

TABLE I 
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TABLE II: Median particle diameter (d), effective anisotropy constant K e ff and energy barrier 
Ea/ks obtained from thermal dependence of coercive field (Hc(T)); Langevin fits of the SPM 
magnetization (Mspm)', thermal dependence of the hyperfine fields (eq|BJ). 

(d) {nm) a K eff (10 4 erg/cm 3 ) b E a /k B (K) 





END1 


ENDS 


END1 


ENDS 


END1 


ENDS 


Hc(T) 


8.6 


9.0 


20.7 


23.5 


449 


509 


MsPM 


7.7 c 


7.9 C 






323 a 


349 a 


Moss 




9.1 




24.3 




526 


Xac(f,T) d 


9.9 


10.3 


32.1 


35.1 


694 


763 



°Values calculated using K^f = 18.7 x 10 4 erg/cm 3 . 



Values calculated using {d) TEM = 8.3 nm. 



c Calculatec 
d From Ref. 



ing the experimental Ms = 7.2 — 7.3 emu/g values at T = 250K. 



TABLE II 
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TABLE III: Hyperfine parameters at different applied fields B app . B e ff is the effective field, It is the total intensity of each six-line 
subespectrum, A is the ratio of intensities of lines 2 and 1 (A2/A1), and (3 is the canting angle of each sublattice. 



B app (kG) 









10 




30 




50 




80 




120 


Sites 


A 


B 


A 


B 


A 


B 


A 


B 


A 


B 


A 


B 


Beff (kG) 


526(1) 


508(1) 


530(1) 


505(1) 


540(1) 


496(1) 


565(1) 


479(1) 


592(1) 


449(1) 


632(1) 


411(1) 


IS (mm/s) 


0.18(1) 


0.24(1) 


0.17(1) 


0.24(1) 


0.15(1) 


0.23(1) 


0.18(1) 


0.26(1) 


0.13(1) 


0.25(1) 


0.13(1) 


0.25(1) 


QS (mm/s) 


0.03(2) 


-0.02(2) 


0.03(2) 


-0.03(2) 


0.00(2) 


-0.04(1) 


0.01(1) 


-0.04(2) 


0.02(1) 


-0.01(2) 


-0.01(1) 


-0.01(1) 


r (mm/s) 


0.63(3) 


0.86(4) 


0.45(4) 


0.79(4) 


0.69(4) 


0.86(3) 


0.62(4) 


0.83(4) 


0.70(4) 


0.82(3) 


0.75(4) 


0.84(4) 


It (%) 


35(2) 


65(2) 


34(3) 


66(2) 


35(3) 


65(2) 


38(2) 


62(3) 


60(2) 


40(3) 


61(2) 


39(3) 


A 


0.70(3) 


0.68(3) 


0.04(3) 


0.23(3) 


0.00(4) 


0.16(4) 


0.02(3) 


0.13(3) 


0.00(3) 


0.09(3) 


0.00(3) 


0.10(3) 


(3 a (degrees) 






0(4) 


34(4) 


0(4) 


31(4) 


0(4) 


31(4) 


0(4) 


22(4) 


0(4) 


22(4) 



Calculated from eqEl 



FIGURE CAPTIONS 



Figure V. Powder x-ray-diffraction profile of the lyophilized sample ENDS indexed with 
the (hkl) reflections of the cubic Fe3C>4 phase. The arrows correspond to organic materials. 
Upper left panel: selected-area from TEM images of Fe3C>4 nanoparticles where some 
surface faceting is observable. Upper right panel: histogram of the particle size populations 
as observed from TEM images. The solid line is the best fit using a log-normal distribution 
with (d) TEM = 8.3 ± 0.8 nm and cr^ = 0.6. 

Figure 2: Magnetization curve of END1 sample after field-cooling {Hpc = 70 kOe) from 
300 K, and measured for increasing T values at the same field. The arrow indicates the 
melting point of the frozen suspension. 

Figure 3: Fit of M(H, 250 K) curve for sample END1 using equation ^ The values of ~p 
and <r M obtained for the best fit are 290/i# and 1.68, respectively. 

Figure 4: Temperature dependence of the coercive field Hc(T) extracted from M(T, H) 
cycles. The solid lines represents the linear fits using eq. El for samples ENDS and END1. 
The inset shows typical magnetization curves M(H, T) for lyophilized sample ENDS taken 
at different temperatures. Note that M(H) approaches to saturation for fields H > 10 kOe. 

Figure 5: Mossbauer spectra of sample ENDS at 4.2 K at different fields B app applied 
along the 7-ray direction. Solid lines are the fitted experimental spectra (open circles), and 
dashed lines correspond to each component with the hyperfine parameters shown in Table ITTT1 

Figure 6: T B (hif, H) dependence extracted from T, H) data taken with H = 10 
and 150 Oe. The solid lines represent the best fit with eq. EJ 

Figure 7: The dependence Tb{H) from ZFC curves of samples ENDS and END1 (I s * 
and 2 nd run) at different fields. 
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